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INTRODUCTION 


An important consideration in the development of modern helicopters is the 
vibratory response of the main rotor blade. One way to minimize vibration levels is 
to ensure that natural frequencies of the spinnning main rotor blade are well 
removed from integer multiples of the rotor speed. This report demonstrates a 
technique for dynamically tuning a finite -element model of a rotor blade to 
accomplish that end. 

Rotor blades are an ideal subject for this type of analysis because a good 
structural representation can be achieved with a single string of beam elements and 
relatively few degrees of freedom. This means that the numerous system stiffness 
and mass matrices required can be formed with relatively low central processor 
time. The technique is valid, however, for larger and more complex models. 

Because the tuning process involves the independent redistribution of mass and 
stiffness, it is especially applicable to composite blade designs in which mass and 
stiffness can be controlled independently by fiber orientation and the use of 
nonstructural mass. 

In the following sections, a brief overview is given of the general purpose 
finite element system known as Engineering Analysis Language (EAL, ref. 1) which was 
used in this work. A description of the EAL System Modification (SM) processor is 
then given along with an explanation of special algorithms developed to be used in 
conjunction with SM. Finally, this technique is demonstrated by dynamically tuning 
a model of an advanced composite rotor blade. 

This work was accomplished in support of the Interdisciplinary Research Office 
of NASA Langley Research Center and the objectives were threefold. The first was to 
establish a technique for tuning the natural frequencies of a spinning rotor blade. 
The second was to demonstrate the usefulness of the EAL SM processor and to be able 
to perform sensitivity and modification operations without dependence on additional 
software. The final objective was to provide guidelines on advanced use of the SM 
processor, i.e., use beyond the scope of currently available documentation. 


ENGINEERING ANALYSIS LANGUAGE (EAL) 


EAL is a general purpose finite element system produced by Engineering 
Information Systems, Inc. It evolved from an earlier finite element program known 
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as SPAR (ref. 2). In its present form, EAL consists of an Executive Control System 
(ECS) in which the user can execute work flow logic, looping, branching and data 
storage; and processors (similar to subroutines) which actually perform structural 
and utility computations. Data input or computations result in data sets which are 
stored in binary data bases or libraries which can be saved and referred to 
indefinitely. The user communicates with and uses these features with input known 
as runstreams. 

Reference 1 is the current EAL reference manual, however, the older SPAR 
reference manual (ref. 2) must be used for the SM processor. EAL version 209 was 
used in this work. 

EAL SYSTEM MODIFICATION (SM) FOR FREQUENCY MODIFICATION 

The approach in modifying frequencies is to first specify a set of target 
(required) eigenvalues corresponding to natural frequencies of the original model. 
Parameters to be changed must be identified along with limits on acceptable changes. 
Sensitivities of the eigenvalues to parameter changes must then be calculated. To 
determine the actual structural changes, the statistical method described in 
reference 3 is used. 

SM operates in 4 phases as described below. The notation used here is 
generally consistent with the SM description contained in reference 2. 

Phase 1 : The differences (Ay) between the eigenvalue targets (Xj) and current 

eigenvalues (X) are calculated. That is: 

ax=x t -x (1) 

Phase 2 : The purpose of phase 2 is to approximate the sensitivities of eigenvalues 

(radians^/sec^) to specified changes in structural parameters which affect 
stiffness and/or mass. These specified changes are known as unit parameters. 

System stiffness change (aK) and mass change (aM) matrices are formed for each unit 
parameter. 

Because the original model eigenvalue solution is based on equation 2 below, 
where Xj is the ith eigenvalue and M, K and Y-j are the system mass, system 
stiffness and the ith mode shape, respectively, then the modified system can be 
described by equation 3. 

x i MY i -KY i =0 (2) 

(X 1 +AX i ) ( M+AM) - (K+AK ) ( Y i +A Y i )=0 ( 3 ) 

With some simplifying assumptions (i.e. changes in mode shapes and products of 
the changes (a's) are very small), a reasonable approximation of eigenvalue 
sensitivity is expressed by equation (4). 

AX i =Y. T AKY.-X i Y T AMY i (4) 

The AK, AM and AAj are therefore the results of phase 2 which is computation- 
ally the most costly phase because the system mass and stiffness matrices must be 
formed for each unit parameter. Computations in the other 3 phases are trivial in 
terms of central processor time. 
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Equation 4 is valid only for a nonspinning structure and must be augmented for 

a spinning structure as described later. 

Phase 3 : The actual structural changes needed to realize the targeted eigenvalues 

are estimated based on equation (5) below which is an adaptation of the work 

presented in reference 3.' 

4P| =[S rr ] [NT] T f[NTl [SJ [ NT] T + [S A ‘[N] ]sx] (5) 

where: 

aP is a set of multipliers which reflects the total estimated structural 
modifications needed in terms of corresponding unit parameters. 

S rr is the covariance or weighting matrix. The diagonal terms, each 

corresponding to the unit parameters in sequence, allow for the relative 
weighting of those parameters. In this application, values are set at unity 
and reset in later iterations if the parameter change limits are being 
exceeded. 

N is a matrix containing reciprocals of the current eigenvalues ( 1 /x -§ ) . 

T is the sensitivity matrix consisting of (ax^'s) with the rows corresponding 

to the number of targets and the columns to the number of unit parameters. 

S ee is the target tolerance matrix associated with acceptable variances of the 

resulting eigenvalues from the targets. 


AX is as described in equation (1). 

The purpose of using this method is to achieve the targeted eigenvalues with 
minimum change to the structure. S rr can be used to influence how much a 
particular unit parameter is changed. For example, a unit parameter which can be 
changed with small penalty or is not likely to exceed the prescribed change limits 
may be assigned an S rr value of 1.0, whereas, a unit parameter which should be 
changed as little as possible may be assigned a value of 0.1. $ ee values normally 
range from 0.0 (when a more exact attainment of the targeted eigenvalue is being 
sought) to 0.1 (when only an approximate result is needed). As described in 
reference 4, S ee values of 0.001 when most $ rr values are 1.0 normally provide 
satisfactory results. 

Phase 4 : Each term of the aP matrix is compared to the parameter change limits data 

set (described below). If any of the limits are exceeded, a aPX matrix is formed 
where the smaller terms (from a P or limits) are used. aPX (aP if no limits were 
exceeded) is then used to actually change the structural parameter data sets of the 
finite-element model. 

To test the results after the completion of phase 4, new mass and stiffness 
matrices must be formed, and the original process of computing mode shapes and 
frequencies is repeated. Normally, two or three iterations are sufficient to 
achieve the desired results if reasonable targets, unit parameters and change 
limits were selected. A complete iteration is the execution of phases 1-4 and 
testing of the results by calculating frequencies of the modified structure. 
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Prior to executing SM, EAL data sets must be established defining the targets, 

parameters, change limits, weighting and target tolerances. The EAL data set names 

for these inputs are given below followed by brief descriptions. 

TVAL - Target (desired) eigenvalues (radian^/sec^) preceded by mode sequence 
numbers. 

PARA - Each PARA data set is a group of changes (incremental element parameter or 
rigid mass) expressed as a fraction of the existing value. Each data set 
is then considered a unit parameter in SM computations. 

SEE - Target tolerance matrix (S ee ). 

SRR - Covariance or weighting matrix (S rr ). 

DPLI - Parameter change limits (minimums and maximums) expressed as multiples (+ or 
-) of unit parameters defined in the PARA data sets. 


AUGMENTATION TO THE SM PROCESSOR 


In this application, it was necessary to develop three algorithms to augment 
the SM processor. These were implemented in the EAL Arithmetic Utility System (AUS) 
processor. The first was to add the centrifugal stiffening effect of the mass 
change (aM) matrices to the sensitivity matrix. The second was to revise the 
weighting matrix (S rr ) when the original values resulted in too many values of the 
change limits data set (DPLI) being violated by the aP matrix, thus causing 
structural changes which were inadequate in achieving targeted results. The third 
was to update the change limits after a complete iteration so that in the next 
iteration, the change limits data set (DPLI), which is based on a fraction of the 
current structural data set values, expresses the same engineering limits in terms 
of mass or stiffness originally intended. 

To correct the sensitivity matrix, an additional system stiffness matrix must 
be formed for each nonzero aM matrix formed in phase 2. This matrix [aKC] reflects 
the centrifugal effect of the spinning aM and is formed using the AUS SPIN command 
to calculate a centrifugal force matrix and the elastic and centripetal 
contributions to stiffness. The Static Solution (SSOL) processor is used to 
calculate deflections due to the centrifugal force. The resulting stresses are 
embedded in the element state data sets by the GSF processor. Geometric 
stiffness changes are then calculted using the KG processor. The elastic, 
centripetal and geometric stiffness contributions are then summed to form [aKC] 
which is used to finalize the sensitivities as follows: 

AX iT0TAL =X i +Y i t AKC ] Y { ^ 

where Xi is given by equation 4. 

The weighting matrix (S rr ) is revised when limits (DPLI) are violated by 
(aP). This is accomplished simply by multiplying each term of the (S rr ) matrix by 
the ratios of corresponding terms of the aPX and aP matrices. That is, 
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(7) 


[ Srr ] nevi'M 


OLD 


fiPX^ 

6p;- 


6PX, 

jp.- 


which has the effect of reducing those S rr terms corresponding to unit parameters 
which are tending to be changed beyond their allowable limits in phase 3. This 
process is repeated until the resulting aPX matrix resulting from phase 4 does not, 
in the judgement of the user, differ too greatly from the aP matrix. If this cannot 
be achieved, the targets may be unachievable based on the selected parameters and 
change limits. 

The updating of the change limits (DPLI) for the subsequent iteration is 
achieved by the following process which updates each term of the DPLI matrix to 
retain the original engineering value. 


[apli] 


NEW 


L 10LD" aPX 1 

'I+f^APX'" 


U 10LD -Af>X l 

"I+?Japx"‘ 


L 20LD“ aPX 2 

“l+f‘APX"‘ 


•• (Etc for each parameter) 


( 8 ) 


where: 


Ap LI NEW = New parameter change limits data set. 
l 10LD ,l 20LD = Old lower limits for parameters 1 and 2. 

^10LD = Old upper limit for parameter 1. 

aPXi, aPX 2 = The final changes for parameters 1 and 2 produced in SM phases 3 
and 4. 

fl. ^2 = The fraction used in defining a unit change for parameters 1 and 2 in 
the PARA data sets. For this process to work, the fraction must be 
uniform within a given PARA data set. 


DEMONSTRATION 


The finite element model (see figure 1) used in this report is based on a 
preliminary design of an advanced composite main rotor blade developed by Mark W. 
Nixon of the U.S. Army Aerostructures Research Group at Langley Research Center. 
Table I gives the mass and stiffness properties of the baseline model which resulted 
from a composite analysis program also developed by Mr. Nixon. Table II provides 
the constraints or parameter changes which cannot be exceeded during the tuning 
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process. These constraints are based on the designer's estimate of what changes can 
be reasonably made without sacrificing the structural integrity or performance of 
the rotor blade. 

Additional constraints on the problem were that bending stiffness, if modified, 
must be changed uniformly over large segments of the blade. The minimum allowable 
mass moment of inertia about the hub was 19000 lb- in- sec 2 for autorotation 
capability. 

The objective of the tuning process was to minimize resonances caused when 
flexible mode frequencies were too close to integer multiples of the rotor speed up 
to eight per revolution (8P). The main rotor speed was 263 RPM (4.3833 HZ) and a 
criterion of at least .2P separation was used. Table III lists the unacceptable 
frequency ranges along with the natural frequencies of the original model and those 
of the modified model following the first and second tuning iterations. 

The overall process which was conducted interactively is depicted in figure 2. 
Figure 3 contains the actual EAL runstreams used in the process. The runstreams in 
combination with this paper and the references should provide adequate guidelines 
for a new SM user. 

Modes 1 and 2 are the flatwise and edgewise rigid body modes, and due to the 
physics of a spinning rotor blade, cannot be significantly altered. Modes 3 through 
7 were therefore targeted for modification. Due to blade twist, modes 3, 5, 6 and 7 
are combined flatwise/edgewise bending modes whereas mode 4 is predominantly 
torsion. It appeared reasonable to drive all of the bending mode frequencies to 
approximately . 25P below the nearest P multiple while allowing the torsion mode to 
remain close to its original frequency. A study of the sensitivities indicated that 
to drive frequencies in opposite directions would have required unacceptably large 
changes in certain parameters. The selected target frequencies are listed in Tables 
III and IV. Table IV also lists all of the SM inputs. 

Results of two complete iterations are summarized in Table III and figure 4. 
Figure 4 shows the ratio of calculated to target frequencies plotted against the 
iteration number ("0" iteration being the original model). A ratio of 1.0 would 
indicate complete convergence with the target value. The first iteration did not 
move all of the frequencies to acceptable ranges (Table III) but did move all of 
them towards the targets as shown in figure 4. The second iteration produced 
frequencies out of the unacceptable ranges and very close to the targeted 
frequencies. The total weight of the blade increased from 250.54 lb to 265.30 lb 
and the mass moment of inertia about the hub increased from 19007.1 to 19780.7 
lb-in-sec 2 . Table V summarizes the final structural properties of the modified 
rotor blade model. 


CONCLUDING REMARKS 

A sensitivity technique useful in minimizing vibrations associated with 
helicopter rotor blades has been demonstrated. This and similar techniques can be 
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effective in achieving desired performance with minimum change to the basic 
structure. This is especially true for spinning structures because centrifugal 
stiffening complicates the intuitive process of changing mass and stiffness to tune 
natural frequencies. 

An advantage of the process described in this report is that the modification 
capability is built into the structural analysis program. This eliminates the need 
for data transfer and development or use of external software. 

The EAL System Modification processor has applications beyond that for which it 
was originally produced and documented, as demonstrated here for a spinning 
structure. As long as the equations for calculating appropriate sensitivities are 
known, structural modification can be computed to achieve any targeted response such 
as mode shapes, static deflections, stress and bending moments and loads due to 
dynamic loads. 
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TABLE I. - MODEL PROPERTIES 


J oi nt 
No. 

J oi nt 
Location 
(z.in) 

Lumped 

Mass 

(lb) 

Lumped 
inertia 
about 
z-axis 
( 1 b- 1 n- sec 2 ) 

1 

0 

0 

.2415 

2 

16.1 

0 

.2700 

3 

18.0 

0 

.0585 

4 

20.0 

0 

.1230 

5 

26.2 

0 

.1860 

6 

32.4 

0 

.1860 

7 

38.6 

0 

.1905 

8 

45.1 

0.32 

.1551 

9 

51.5 

0.645 

.1161 

10 

58.0 

0.645 

.1161 

11 

64.4 

1.93 

.3474 

12 

96.6 

3.22 

.5796 

13 

128.8 

3.22 

.5796 

14 

161.0 

3.22 

.5796 

15 

193.2 

2.415 

.4347 

16 

209.3 

1.61 

.2898 

17 

225.4 

1.61 

.2898 

18 

241.5 

1.61 

.31395 

19 

257.6 

1.61 

.36225 

20 

273.7 

1.61 

.3864 

21 

289.8 

1.125 

.2700 

22 

296.2 

0.645 

.1548 

23 

302.7 

0.645 

.1548 

24 

309.1 

0.32 

.2068 

25 

315.6 

0.35 

.2580 

26 

322.0 

0.35 

.1280 






Edgewise 

Flatwise 

Twist 


Cross 

Torsional 





Stiffness 

Stiffness 

Angle 

Distributed 

Sectional 

Stiffness 

Beam 

Joi nts 

EIll* 

El 22* 

LE Down 

Weight 

Area 

GJ 

Section 

Spanned 

(LBF-in 2 ) 

(LBF-in 2 ) 

(DEGR) 

(Lbs/in) 

In 2 

(LBF-in 2 ) 

1 

1 

_ 

2 

900.0 

900.0 

26.0 

2.29 

44.44 

100.00 

2 

2 

- 

3 

.0001 

.0001 

25.34 

2.29 

44.44 

87.50 

3 

3 

- 

4 

.0001 

.0001 

25.26 

2.29 

44.44 

87.50 

4 

4 

- 

5 

580.0 

360.0 

25.17 

2.20 

44.44 

75.50 

5 

5 

- 

6 

580.0 

360.0 

24.92 

2.20 

44.44 

75.50 

6 

6 

- 

7 

580.0 

360.0 

24.66 

2.20 

44.44 

75.50 

7 

7 

- 

8 

580.0 

298.0 

24.41 

2.60 

44.44 

60.00 

8 

8 

- 

9 

1260.0 

25.89 

24.14 

0.35 

78.84 

17.125 

9 

9 

- 

10 

1260.0 

25.89 

23.87 

0.35 

78.84 

17.125 

10 

10 

- 

11 

1260.0 

25.89 

23.60 

0.35 

78.84 

17.125 

11 

11 

- 

12 

1260.0 

25.89 

23.34 

0.35 

78.84 

17.125 

12 

12 

- 

13 

1260.0 

25.89 

22.01 

0.35 

78.84 

17.125 

13 

13 

- 

14 

1260.0 

25.89 

20.68 

0.35 

78.84 

17.125 

14 

14 

- 

15 

1260.0 

25.89 

19.35 

0.35 

78.84 

17.125 

15 

15 

- 

16 

1260.0 

25.89 

18.02 

0.35 

78.84 

17.125 

16 

16 

- 

17 

1260.0 

25.89 

17.34 

0.35 

78.84 

17.125 

17 

17 

- 

18 

1260.0 

25.89 

16.69 

0.35 

78.84 

17.125 

18 

18 

- 

19 

1260.0 

25.89 

16.03 

0.35 

78.84 

17.125 

19 

19 

- 

20 

1260.0 

25.89 

15.67 

0.35 

78.84 

17.125 

20 

20 

- 

21 

1260.0 

25.89 

14.70 

0.35 

78.84 

17.125 

21 

21 

- 

22 

1260.0 

25.89 

14.03 

0.35 

78.84 

17.125 

22 

22 

- 

23 

1260.0 

25.89 

13.77 

0.35 

78.84 

17.125 

23 

23 

- 

24 

1260.0 

25.89 

13.51 

0.35 

78.84 

17.125 

24 

24 

- 

25 

580.0 

24.0 

13.23 

0.90 

44.44 

60.00 

25 

25 

- 

26 

580.0 

24.0 

12.77 

3.3706 

44.44 

60.00 

Total 

mass moment of 

inertia about x 

axis (hub): 19007.1 

lb in sec 2 . 

Total weight: 

250.54 lb. 



* Stiffness paramaters are with respect to a local reference frame which is rotated the amount of the twist angle 
from the global frame shown on figure 1. 
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TABLE II. - CHANGE LIMITS ORIGINAL PAGE JS 

OF POOR QUALITY 



Lumped 



Edgewise 


Joint 

Mass 


Section 

Stiffness Flatwise 

No. 

% 


No. 

% 

Stiffness 

1 

0 


1 

0 

0 


2 

0 


2 

0 

0 


3 

0 


3 

0 

0 


4 

0 


4 

0 

0 


5 

0 


5 

0 

0 


6 

0 


6 

0 

0 


7 

0 


7 

0 

0 


8 

0 


8 

+ 20 " 

+10 n 


9 

-50 +100 


9 

+20 

+10 


10 

-50 +100 


10 

+20 

+10 


11 

-50 +100 


11 

+20 

no 


12 

-50 +100 


12 

+20 

+10 

13 

-50 +100 


13 

+20 

no 


14 

-50 +100 


14 

+20 

+10 


15 

-50 +100 


15 

+20 

>* +10 J 


16 

-50 +100 


16 

+20 

+10 


17 

-50 +100 


17 

+20 

no 


18 

-50 +100 


18 

+20 

no 


19 

-50 +100 


19 

+20 

no 

* 

► 

20 

-50 +100 


20 

+20 

no 

21 

-50 +100 


21 

+20 

no 


22 

-50 +100 


22 

+20 

no 


23 

-50 +100 


23 

+20J 

+10 > 


24 

-50 +100 


24 

0 

0 


25 

+100 


25 

0 

0 


26 

+100 






Mi nimum 
* Items 

Allowable Mass Moment of Inertia About X Axis: 19000 lb ■ 

in brackets must be changed uniformly as a group. 


TABLE III. - MODEL NATURAL FREQUENCIES 
COMPARED TO UNACCEPTABLE RANGES 


MULTIPLE, M UNACCEPTABLE RANGES, HZ (MP±.2P) 


1 

3.507 - 5.260 

2 

7.890 - 9.643 

3 

12.273 - 14.027 

4 

16.657 - 18.410 

5 

21.040 - 22.793 

6 

25.423 - 27.177 

7 

29.807 - 31.560 

8 

34.190 - 35.943 


WHERE P=263rpm OR 4.3833HZ 


FREQUENCIES 


MODE 

TARGET 

ORIGINAL 

ITER 1 

ITER 2 

3 

12.054 

12.488* 

12.067 

12.030 

4 

16.0896 

16.090 

16.010 

15.988 

5 

20.8208 

22.460* 

21.062* 

20.928 

6 

23.0125 

25.056 

23.158 

22.949 

7 

33.9708 

36.368 

34.326* 

34.114 


*IN UNACCEPTABLE RANGE 


255 


TABLE IV 


SYSTEM MODIFICATION INPUT DATA 


TARGET FREQUENCIES 
(TVAL SM) 1 


SEQUENCE NO. 

MODE NO. 

EIGENVALUE 



1 

3 

5736.3299 

(12.05 

HZ) 

2 

4 

10219.9920 

(16.09 

HZ) 

3 

5 

17114.1744 

(20.82 

HZ) 

4 

6 

20906.7892 

(23.01 

HZ) 

5 

7 

45558.7855 

(33.97 

HZ) 


UNIT PARAMETERS 
(PARA SM n) 



DATA 

LINE 2 


COLUMN 3 

n 

TYPE 

NO. 

FRACTION 

NO. 

1 

RIGID MASS 

25 

.1 

1 TO 3 


RIGID MASS 

26 

.1 

l TO 3 

2 

RIGID MASS 

9 

.1 

1 TO 3 


RIGID MASS 

10 

.1 

1 TO 3 

3 

RIGID MASS 

11 

.1 

1 TO 3 

4 

RIGID MASS 

12 

.1 

1 TO 3 

5 

RIGID MASS 

13 

.1 

1 TO 3 

6 

RIGID MASS 

14 

.1 

1 TO 3 

7 

RIGID MASS 

15 

.1 

1 TO 3 

8 

RIGID MASS 

16 

.1 

1 TO 3 

9 

RIGID MASS 

17 

.1 

1 TO 3 

10 

RIGID MASS 

18 

.1 

1 TO 3 

11 

RIGID MASS 

19 

.1 

1 TO 3 

12 

RIGID MASS 

29 

.1 

1 TO 3 

13 

RIGID MASS 

21 

.1 

1 TO 3 

14 

RIGID MASS 

22 

.1 

1 TO 3 


RIGID MASS 

23 

.1 

1 TO 3 

15 

RIGID MASS 

24 

.1 

1 TO 3 

16 

EDGEWISE STIFFNESS 





< EI 11> 

8 TO 23 

.1 

4 

17 

FLATWISE STIFFNESS 





(EI 22 ) 

16 TO 23 

.1 

6 

18 

FLATWISE STIFFNESS 





(EI 22 ) 

8 TO 15 

.1 

6 


TARGET 

NO. 

TARGET TOLERANCES 
(SEE SM) 

MODE 

NO. 

TOLERANCE 

1 

2 

3 

4 

5 

UNIT 

PARAMETER NO. 

3 

4 

5 

6 
7 

INITIAL COVARIANCE 5 
AND CHANGE LIMITS 
(SRR SM AND DPLI SM) 

COVARIANCE 

.001 

.1“ 

.001 

.001 

.001 

LIMITS 6 

1 

1 

-10, +10 

2-15 

1 

-5, +10 

16 

1 

-2, +2 

17 

1 

-1,+1 

18 

I 

-1,+1 


NOTES: These data correspond to the input in the EAl runstream in figure 3b. 

1. Names in parentheses are EAL data set names. 

2. Line number of structural data set corresponds to joint for rigid 
masses and beam segment number for stiffnesses. 

3. The unit parameter is a set of numbers computed from multiplying 
the fraction times the structural values in the indicated lines 
and columns. 

4. A tolerance value of 0.1 rather than 0.001 indicates that it is 
less critical for the final frequency to be very close to the 
target value. 

5. These values were modified in the iteration process. 

6. Limits of -5 to plus 10 means that the structural parameter cannot 
be reduced by more than 5 x (FRACTION) x (EXISTING VALUE) nor 
increased more than 10 x (FRACTION) x (EXISTING VALUE). 


TABLE V. 

- FINAL 

MODIFIED 

STRUCTURAL 

PROPERTIES 


LUMPED 


EDGEWISE 

FLATWISE 

JOINT 

MASS 

BEAM 

STIFFNESS 

STIFFNESS 

NO. 

(lb) 

SECTION 

EIn , 

ei 22. 




(LBF in 2 ) 

LBF in 2 ) 

9 

1.284 

8 

1046.96 

23.311 

10 

1.284 

9 

1046.96 

23.311 

a 

3.785 

10 

1046.96 

23.311 

12 

6.440 

11 

1046.96 

23.311 

13 

6.286 ! 

12 

1046.96 

23.311 

14 

6.401 

13 

1046.96 

23.311 

15 

1.883 

14 

1046.96 

23.311 

16 

1.025 

15 

1046.96 

23.32 

17 

1.881 

16 

1046.96 

23.32 

18 

3.190 

17 

1046.96 

23.32 

19 

3.220 

18 

1046.96 

23.32 

20 

3.201 

19 

1046.96 

23.32 

21 

0.5713 

20 

1046.96 

23.32 

21 

0.3225 

21 

1046.96 

23.32 

23 

0.3225 

22 

1046.96 

23.32 

24 

0.1815 

23 

1046.96 

23.32 

25 

0.1275 




26 

0.1275 





TOTAL MASS MOMENT OF INERTIA ABOUT X AXIS (HUB): 19780 lb-in sec 2 . 


TOTAL WEIGHT: 265.30 lb 

NOTE: All other properties unchanged from Table I. 
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FINITE ELEMENT MODEL 




200 IN 300 IN 


Z 


CROSS SECTION 



Figure 1.- Rotor blade model. 



Figure 2.- Procedure for tuning frequencies of a 
spinning structure. 



*XQT TAB 
START 26 

JLOC 

1 0. 0. 00. 

2 0. 0. 16.1 

3 0. 0. 18.0 

4 0. 0. 20.0 

5 0. 0. 26.20 

6 0 . 0 . 32.40 

? 0. 0. 38.60 

8 0. 0. 45.10 

9 0. 0. SI. 50 

10 0. 0. 58.00 

11 0. 0. 64.40 

12 0 . 0 . 96.60 

13 0. 0. 128.86 

14 0. 0. 161.06 

15 0. 0. 193.20 

16 0. 0. 209.36 

17 0. 0. 225.40 

18 0. 0. 241.56 

19 0. 0. 257.60 

20 0. 0. 273.70 

21 0. 0. 289.80 

22 0. 0. 296.29 

23 9. 0. 302.70 

24 0. 0. 309.10 

25 0. 0. 315.60 

26 0. 0. 322.00 

CON- 1 

ZERO 1 2 3 4 5 6U 
fMTC 

t 1.46 .2S 0. BG-.4+6 

BA 

00 I UN II A1 12 A2 A F 
CIUN 1 960.0 0 . 900.0 0 . 44.44 100.00 

CIUN 2 .0001 0. .0001 0. 44.44 87.500 

GIUN 3 .0001 0. .0001 0. 44.44 87.500 

CIUN 4 580.0 0. 360.0 0. 44.44 75.000 

CIUN 5 580.0 0. 360.0 0. 44.44 75.000 

GIUN 6 580.0 0. 360.0 0. 44.44 75.000 

CIUN 7 580.0 0. 298.0 0. 44.44 60.000 

GIUN 8 1260. 0. 25.89 0. 78.84 17.125 

GIUN 9 1260. 0. 25.89 0. 78.84 17.125 

CIUN 10 1269. 0. 25.89 9. 78.84 17.125 

CIUN 11 1269. 0. 25.89 0. 78.84 17.125 

CIUN 12 1260. 0. 25.89 0. 78.84 17.12S 

CIUN 13 1260. 0. 25.89 0. 78.84 17.125 

GIUN 14 1269. 0. 25.89 0. 78.84 17.125 

GIUN 15 1260. 9. 25.89 0. 78.84 17.125 

CIUN 16 1260. 0. 25.89 9. 78.84 17.125 

CIUN 17 1260. 9. 25.89 0. 78.84 17.125 

GIUN 18 1269. 9. 25.89 0. 78.84 17.125 

CIUN 19 1269. 0. 25.89 0. 78.84 17.125 

GIUN 20 1269. 9. 25.89 0. 78.84 17.125 

CIUN 21 1260. 0. 25.89 0. 78.84 17.125 

CIUN 22 1260. 9. 25.89 9. 78.84 17.125 


CIUN 23 1260. e. 25.89 0. 
CIUN 24 580.0 0. 24.00 0. 
CIUN 25 589.9 0. 24.00 0. 
HREF 


1 

1 

1 

-1 

89879 

2 

1 

1 

- 1 

90378 

3 

1 

1 


90438 

4 

1 

1 

-1 

90S05 

5 

1 

1 

-j 

90690 

6 

1 

1 

-1 

99880 

7 

l 

1 

-1 

91061 

8 

1 

1 

- 1 

91255 

9 
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1 

-1 

91447 

10 

1 

1 

-1 

91636 

11 

1 

1 

-I 

91817 

12 

1 

1 

-1 

.92712 

13 

1 

1 

-1 

93557 

14 

1 

1 

-1 

.94351 

15 

1 

1 

”1 

95995 

16 

1 

1 

-I 

.95455 

17 

1 

1 

-1 

.95787 

18 

1 

1 

-J 

.96112 

19 

1 

1 

-1 

.96281 

20 

1 

1 

-1 

.96727 

21 

1 

1 

-1 

.97817 

22 

1 

1 

-I 

.97126 

23 

1 

1 

-1 

.97234 

24 

1 

1 

-1 

.97346 

25 

1 

1 

-1 

•97527 

RNASS 




cn 

0O259067 .0025 

8 


0 

32 


9 


0 

645 



10 0.645 

11 1.93 

12 3.22 

13 3.22 

14 3.22 

15 2.415 

16 1.61 

17 1.61 

18 1.61 

19 1.61 

29 1.61 

21 1.125 

22 9.645 

23 9.645 

24 9.32 

25 .35 

26 .35 

NSU0 NON STRUCT UT DISTR 

1 2.29 

2 2.29 

3 2.29 

4 2.20 

5 2.20 

6 2.29 


78.84 17.125 
44.44 69.000 
44.44 60.090 


7 2.69 

8 0.35 

9 0.35 

10 9.35 
It 9.35 

12 0.35 

13 0.35 

14 0.35 

15 0.35 

16 0.35 

17 0.35 

18 0.35 

19 0.35 

20 0.3S 

21 0.35 

22 0.35 

23 0.35 

24 9.90 

25 3.3706 
tXQT ELD 
E21 

CROUP 1' ROTOR BLADE 

NSECT-1 :NNAT"i:NN$U«l 

INC NSECT-1 

INC NNSU-1 

INC NREF-1 

1 2 1 2S 

tXQT E 

RESET G-386. 

tXQT AUS 

TABUNI-6 HJ-26MIN6 


.241500000 
. 270000000 
.585000E-01 
. 123009000 
. 1860OOO00 
. 186000000 
. 190500009 


1-6* J- 8> 

1-6) J* 9< 

I -6i J- 19> 

I-6i J« ill 

1-61 J- 12t 

1-61 J- 131 

I-fil J- Ml 

1-61 J- 15 * 

I-6i J- 161 

1-61 J- 171 

I-6t J- 18i 

l-6i J- 19i 

i -6i j- an 

I-6i J- 21i 

I-6i J- 22> 

I*6i J- 23i 

I-6i J- 241 

I -6 1 J« 25t 

I-6< J- 261 



DEf12-SUNCDEf1 IN6) 

n+Rn • sumRNAS de«2 j 

tXQT DCU 
DISA 1 EQNF 
tXQT EKS 
tXQT TAN 
tXQT K 
SPDP-2 
tXQT AUS 

SPIN! H+RN K 27.5413 9.9 0. 9. 9. 0. S 27.5413 RAD'SEC 

tXQT RSI 

RESET K-KSPN 

tXQT SSOL 

RESET K-KSPN 

tXQT CSF 

RESET ENBCD-1 

tXQT ICC 

SPDP-2 

tXQT AUS 

ICECG-SUmiCSPN ICC) 
tXQT RSI 
RESET IC-KECC 
tXQT EIG 

RESET INIT-11 NREQ-7 N-H+RN IC-ICECG OUTL-1 
tXQT AUS 
RIC-RICIDC 1 ) 

NR-PROD(N*RN RIG) 

GN-XTYD(RIC NR) 

CHUT -UN I ON (386. CfU 
tXQT DCU 
PRINT 1 GN 
PRINT 1 GPHJT 
tXQT EXIT 



Figure 3a.- EAL runstream for calculation of natural frequencies 
of the spinning structure. 
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IXQT AUS 

TABL(NI-2 NJ-SUTUAL SPI 
J • 1 * 3 . 5736.3298 $12,654 HZ 
J*2M. 10219.9929 9 16.0896 
J-3iS. 17114.1744 $ 20.8298 
J-4I6. 20906.7892$ 23.0126 
J'5t7. 45558.7855$ 33.9798Z 


TABLtNI-5 

NJ-2 ) iPARA 

Srt 

liJ-1118. 25. 

.1 

1. 

3 



-2U8. 26. .1 1. 

3. 



TABL(NI-S 

N J-2): PARA 

Srt 

2IJ-1I18. 09. 


1 , 

3 




J-2118. 19. 


1. 

3 

TABL(NI-5 

NJ-l ) IPARA 

sn 

3: J- 1 : 18 . 11. 


1. 

3 

TABKNI-5 

NJ-l ) iPARA 

sn 

4IJ-1I18. 12. 


1. 

3 

TABL<NI-S 

NJ-l ) IPARA 

Srt 

5U-U18. 13. 


1. 

3 

TABKNI-5 

NJ-l) IPARA 

Srt 

6: J-l : 18. 14. 


1. 

3 

TABL(NI-5 

NJ-l ) IPARA 

Srt 

7: J-l : 18. 15. 


1. 

3 

TABL(NI«5 

NJ-l )IPARA 

Srt 

8: J-l : 18. 16. 


1. 

3 

TABL(NI»5 

NJ-l ):PARA 

SH 

09U-1U8. 17. 

, . 

1 1 



9. 13. .1 6. 6. 

9. 14. .1 6. 6. 

9. 15. .1 6. 6. 

TABU N I ■ 1 NJ-5USEE SHU-1 5:0.001 .1 
TABU NX -l NJ-18 > > SRR SHiJ-1 18U. 
TABUNI-2 NJ-18 ) tDPLl SrtiJ-ll-10. 10. 

J-2 151-5.0 19. 
J-16I-2. 2. 

J-171-1. |. 

J-181-1. 1. 

IXQT EXIT 


TABUNI 

•5 NJ- 

1 ) IPARA 

Srt 

101 J-l 1 18. 

TABLCNI 

•5 NJ- 

DIPARA 

Srt 

lltj-l t 18. 

TABLtNI 

•5 NJ- 

1 ) i PARA 

Srt 

12: J-l 1 18. 

TABUNI 

•5 NJ- 

1 ) IPARA 

sn 

13: J-l : 18. 

TABUNI 

-5 NJ- 

2) iPARA 

SH 

14t J-l 1 18. 





J-2U8. 

TABUNI 

•5 NJ- 

1 ) IPARA 

SH 

15: J-l t 18. 

TAIL(NI 

•5 NJ- 

16 ) iPARA Srt 16IJ-1 16 

9. 8. 

.1 4. 

4. 



9. 9. 

.1 4. 

4. 



9. 10. 

.1 4, 

4. 



9. 11. 

.1 4; 

4. 



9. 12. 

.1 4. 

4. 



9. 13. 

.1 4. 

4. 



9. 14. 

.1 4 

. 4. 



9. 15. 

.1 4 

. 4. 



9. 16. 

.1 4. 

4. 



9. 17. 

.1 4. 

4. 



9. 18. 

.1 4 

. 4. 



9. 19. 

.1 4. 

4. 



9. 20. 

.1 4. 

4. 



9. 21. 

.1 4. 

4. 



9. 22. 

.1 4. 

4, 



9. 23. 

.1 4 

. 4. 



TABL(NI 

•5 NJ- 

8) iPARA 

Srt 

17i J- 1 8 

9. 16. 

.1 6. 

6. 



9. 17. 

.1 6. 

6. 



9. 18. 

.1 6 

. 6. 



9. 19. 

.1 6. 

6. 



9. 20. 

.1 6. 

6. 



9. 21. 

.1 6. 

6. 



9. 22. 

.1 6. 

6. 



9. 23. 

.1 6 

. 6. 



TABLtNI 

•5 NJ- 

8) IPARA 

sn 

18* J-l 8 

9. 8. 

.1 6. 

6. 



9. 9. 

.1 6. 

6. 



9. 10. 

.1 6. 

6. 



9. 11. 

.1 6. 

6. 



9. 12. 

.1 6. 

6. 




.1 1. 3. 
.1 1. 3. 
.1 1. 3. 
.1 I. 3. 
.1 1. 3. 
.1 1. 3. 


Figure 3b.- EAL runstream for input to the SM processor. 


IXQT SH $ DEVELOP DH FOR EACH PARAH 
RESET NPARA-18 G-386. OUTL-1 HUUX-2 
OPER 1100 
IXQT RSI 

RESET K »K$LAST KECG 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PAR Art 1 
DEFINE DrtT-DH SH 1 1 iDHA-UNION( DHT ) 

SPIN i DrtA K 27.5413 0. 0. 0. 8. 0. 

IXOT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET EMBED- 1 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 0 DK1 SPAR 36 1 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PARArt 2 

DEFINE DHT-DH SH 2 1 : DMA -UN ION ( DHT ) 

SPimDrtA K 27.5413 0. 0. 0. 0. 0. 

IXQT SSOL . 

RESET K-KSPN 
IXQT GSF 
RESET EMBED -1 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 0 DK2 SPAR 36 2 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PARArt 3 

DEFINE DNT-Drt Srt 3 1 t DrtA -UN I ON (DHT) 

SPIN: DrtA K 27.5413 0. 9. 0. 9. 9. 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET EMBED- 1 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 0 1*3 SPAR 36 3 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PARArt 4 

DEFINE DHT-DH Srt 4 1 :DHA-UNION(DHT ) 

SPIN: DMA K 27.5413 0. 9. 0. 0. 0. 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET EMBED- 1 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 0 DK4 SPAR 36 4 

IXQT AUSf DEUEL INERTIAL STIFFNESS FOR PARAH 5 

DEFINE DHT-DH SH 5 1 : DrtA- UN ION (DHT > 

SPIN: DMA K 27.5413 9. 9. 9. 0. 0. 

IXQT SSOL . 

RESET K-KSPN 
IXQT GSF 
RESET EMBED- 1 
IXQT KG 
IXQT XU 

CHAN 1 KG SPAR 36 9 DK5 SPAR 36 5 

IXQT AUSI DEUEL INERTIAL STIFFNESS FOR PARArt £ 

DEFINE DrtT-DH Srt 6 1 iDHA-UNIONCDHT ) 


SPIN! DrtA K 27.5413 9. 9. 9. 0. 0. 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ErtBED-1 
IXOT KG 
IXQT XU 

CHAN 1 KG SPAR 36 0 DK6 SPAR 36 6 

IXQT AUS0 DEUEL INERTIAL STIFFNESS FOR PAR** 7 

DEFINE DHT-DH SH 7 1 :DMA-UNION( DrtT ) 

SPIN:DHA K 27.5413 9 . 0 . 0 . 0 . 9 . 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ErtBED-1 
IXQT KG 
IXQT XU 

CHAN 1 KG SPAR 36 0 DK7 SPAR 36 7 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PARAH 8 

DEFINE DrtT-Drt Srt 8 l : DMA- UN ION (DrtT ) 

SPIN:DHA K 27.5413 0. 0. 9. 9. 0. 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ErtBED-1 
IXQT KG 
IXQT XU 

CHAN 1 KG SPAR 36 9 DC 8 SPAR 36 8 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PARAH 9 

DEFINE DrtT-Drt Srt 9 1 : DMA- UN ION ( DHT ) 

SPINiDrtA K 27.5413 0. 9. 0. 0. 0. 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ErtBED-1 
IXQT KG 
IXQT XU 

CHAN 1 KG SPAR 36 9 DK9 SPAR 36 9 

IXQT AUS$ DEUEL INERTIAL STIFFNESS FOR PARArt 10 

DEFINE DrtT-Drt Srt 19 l : DMA -UN ION (DHT) 

SPIN* DrtA K 27.5413 9 . 0 . 0 . 0 . 0 . 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ErtBED-1 
IXQT KG 
IXOT XU 

CHAN 1 KG SPAR 36 9 DC 19 SPAR 36 19 

IXQT AUSf DEUEL INERTIAL STIFFNESS FOR PARArt 11 

DEFINE DHT-DH SH 11 1 iDrtA-UNION(DHT) 

SPIN (DMA K 27.5413 0 . 9 . 9 . 9 . 0 . 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ERBED-1 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 0 DCU SPAR 36 11 


Figure 3c.- EAL runstream for developing the sensitivity matrix. 
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tXQT AUS* DEUEL INERTIAL STIFFNESS FOR PARAN 12 
DEFINE DNT-DH SH 12 1 tDHA-UNION< DHT > 

SPINtDNA K 27.5413 0. 0. 0. 0. 0. 

tXOT SSOL 

RESET K-KSPN 

IXQT GSF 

RESET EHBED-1 

IXQT KG 

IXQT DCU 

CHAN 1 ICG SPAR 36 0 DK12 SPAR 36 12 

IXQT AUS0 DEUEL INERTIAL STIFFKSS FOR PARAN 13 

DEFINE DNT-DN Sfl 13 1 tDHA-UNION(DRT) 

SPINtDNA K 27.5413 0. I. 0. 0. 0. 

IXQT SSOL 
RESET K-KSPN 
tXQT GSF 
RESET EHBED-t 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 9 DK13 SPAR 36 13 

IXQT AUS* DEUEL INERTIAL STIFFNESS FOR PARAN 14 

DEFINE DNT-DN SH 14 ltDHA-UNION(DHT) 

SPINtDNA K 27.5413 0. 0. 0. 0. 0. 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ENDED* 1 
IXQT KG 
IXQT DCU 

CHAN 1 KG SPAR 36 0 DK14 SPAR 36 14 


IXQT AUS* DEUEL INERTIAL STIFFNESS FOR PARAN 15 
DEFINE DNT-DH SN 15 1 tDNA-UNION(DNT) 

SPINiDRA K 27.5413 0 . 0 . 0 . 0 . 0 . 

IXQT SSOL 
RESET K-KSPN 
IXQT GSF 
RESET ENDED- 1 
IXQT KG 
IXQT DCU 


CHAN 1 KG SPAR 36 0 DK15 SPAR 36 15 


fXXQT AUS* DEUEL INERTIAL STIFFNESS FOR PARAN 16 
•DEFINE DNT-DH SH 16 1 *DHA-UNlON<DNT ) 

•SPINtDNA K 27.5413 0 . 0 . 0. 0 . 0. 

•IXQT SSOL 
•RESET K-KSPN 
•IXQT GSF 
•RESET ENDED- 1 
•IXQT KG 
•IXQT DCU 

•CHAN 1 KG SPAR 36 0 DK16 SPAR 36 16 

•IXQT AUS* DEUEL INERTIAL STIFFNESS FOR PARAN 17 

•DEFINE DNT-DH SH 17 1 tDHA-UNION(DHT) 

•SPINtDNA K 27.5413 *. 0 . 0. 0 . 0 . 

•IXQT SSOL 
•RESET K-KSPN 
•IXQT GSF 
•RESET ENDED- 1 


•IXQT KG 
•IXQT DCU 

•CHAN 1 KG SPAR 36 0 DK17 SPAR 36 17 
IXQT AUS 

DEFI UN5-UIBR NODE 1137 
VK1 -PROD CDK1 UNS) 

YK2-PR0DCDK2 UN5) 

YK3-PR0DC DK3 UH5 ) 

YK4-PR0DCDK4 UN5 ) 

YK5-PR0D( DK5 UNS) 

YK6-PROD(DK6 UNS) 

YK7-PR0DCDK7 UIS) 

YK8-PR0D( DK8 UNS) 

YK9-PR0D( DK9 URS> 

YK10-PRODCDK10 UNS) 

YK11-PRODCDK11 UNS) 

YK12-PR0DCDK12 UNS) 

YK13-PR0D(DK13 UNS) 

YK14-PR0D(DK14 UNS) 

YK1S-PR0D< DK15 UNS) 

SYK16-PR0D(DK16 UN5 ) 

•YK 17-PROD ( OKI 7 UNS) 

SE01-XTYD( UN5 VK1 ) 

SE02-XTYDCUN5 VK2 ) 

SE03-XTYD(UH5 YK3 ) 

SE04-XTYD(UN5 VK4 ) 

SE06-XTYD(UN5 VK5 ) 

SE06-XTYD( UN5 YK6 ) 

SE07-XTYO(UN5 YK? ) 

SE08-XTYD(UNS VK8 ) 

SE09-XTVD(UN5 YK9 ) 

SE10-XTYD(UN5 YK10) 

SE1 1 -XTVD<UN5 YK11) 

SE12-XTYD(UNS VK12) 

SE13-XTYD(UNS YK13) 

S£14«XTYD( UNS VK14) 

SE15-XTYD(UN5 YK1S) 

*SE16-XTYD(UNS VK16) 

•SE17-XTVD(UN5 YK17) 

TABLCNI-S NJ-1 USE16 
TABUNI-5 N J • 1 ) t S€ 1 7 
TABKNI-S NJ-1 USE 18 

SENK-UNIONCSE01 SE62 SE03 SE04 SE05 SE06 SE07> 
SERB SE09 SE10 SE11 SE12 SE13> 

SE14 SE15 SE16 SE17 SE18) 

SEND - SUN (SENK SENS) 

IXQT DCU 

CHAN 1 SEND AUS 1 1 SENS NATR 0 1 
IXQT AUS 

DEFI DA9-DA BTAB 2 9<D9-UNI0N( BA9 ) 

DEFI RH18-RNAS ITAB 2 lBtR18-UNI0N(RH18 ) 

IXQT SN 

RESET OUTl-2 HUUX-2 NUDP-2 NPARA-18 G-3D6. 

OPER 0011 
IXQT AUS 

BA BTAB 2 9-UNI0N(B9) 

RNAS BTAB 2 18-UNION(R10 ) 

IXQT EXIT 


Figure 3c.- EAL runstream for developing the sensitivy matrix (concluded). 



Figure 3d.- EAL runstream for revision of the covariance (weighting) 
matrix. 
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RATIO CALCULATED/TARGET FREQUENCIES 


tXOT Sfl 

RESET OUTl-2 NUUX-2 NUDP-2 NPARA-18 G-386. 

oper e ® ® 1 

TXQT AUS 

TABUNI-l NJ-18)*SRR SHiJ-1 18U. 

DEFI DPSN 2 DPX REV 
DPST -RTRAN< DPSN) 

TABUNI-2 NJ-18 ) tDPCH 

TRAM ( SOUR- DPST ILIN-1 JLIN-18 DSKIP-1 SBASE-® DBASE-®) 
TRAM ( SOUR -DPST 1UM-1 JUN-18 DSKIP-1 SBASE-B DBASE-1) 
DEFI DPO-DPLI 
DPME-SUH(DPO -1. DPCH ) 

DPFR-SUN(SRR .1 DPSP) 

DPRE-RECIP(DPFR) 

TABUNI-2 NJ-18)tRADP 

TRAN < SOUR -DPRE IUN-1 JlIN-18 DSKIP-1 SBASE-B DBASE-® > 
TRAN ( SOUR -DPRE ItIN-1 JLIN-18 DSKIP-1 SBASE-9 DBASE- 1) 
DPN2*PR0D( RADP DPNE ) 

DPLI SN-UNI0WDPN2) 

tXQT E 

RESET G-386. 

tXQT EKS 

IXOT K 

tXQT AUS 

NNEU-$UH(DEN2 RNAS ) 

NDIF-SUmNNEU -1. fHRN) 

N+RN-SUN< DEN2 RNAS) 

SPINiN+RN K 27.5413 ®. ®. ®. ®. ®. 

tXOT RSI 

RESET K-KSPN 

tXOT SSOL 

RESET K-KSPN 

tXOT GSF 

RESET ENDED- 1 

tXQT KG 

tXQT AUS 

KECG-SUN(KSPM KG) 
tXOT RSI 
RESET K-KECG 
tXOT EIG 

RESET INLIB-1 N-NfrfiN K-KECG 
tXQT AUS 
RIG-RIGIDU ) 

NR-PROD(N+RN RIG) 

GN-XTVDCRIG NR) 

GNUT-UNI0NO86. GN) 

RNUT -UNION! 386 . RNAS) 

tXOT DCU 

PRINT l GN 
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Figure 3e.- EAL runstream for calculation of frequencies of the 
modified structure and recomputation of the change 
limits for the next iteration. 
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Figure 4.- Results of modification 


